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The hybrid Monte Carlo (HMC) algorithm is applied for the Bayesian inference of the 
stochastic volatility (SV) model. We use the HMC algorithm for the Markov chain Monte 
Carlo updates of volatility variables of the SV model. First we compute parameters of the 
SV model by using the artificial financial data and compare the results from the HMC 
algorithm with those from the Metropolis algorithm. We find that the HMC algorithm 
decorrelates the volatility variables faster than the Metropolis algorithm. Second we 
make an empirical study for the time series of the Nikkei 225 stock index by the HMC 
algorithm. We find the similar correlation behavior for the sampled data to the results 
from the artificial financial data and obtain a value close to one {ip ~ 0.977), which 
means that the time series has the strong persistency of the volatility shock. 

Keywords: Hybrid Monte Carlo Algorithm, Stochastic Volatility Model, Markov Chain 
Monte Carlo, Bayesian Inference, Financial Data Analysis 



1. Introduction 

Many empirical studies of financial prices such as stock indexes, exchange rates 
have confirmed that financial time series of price returns shows various interesting 
properties which can not be derived from a simple assumption that the price re- 
turns follow th e ge ometric Brownian motion. Those properties are now classified 
as stylized facta^^. Some examples of the stylized facts are (i) fat-tailed distribu- 
tion of return (ii) volatility clustering (iii) slow decay of the autocorrelation time 
of the absolute returns. The true dynamics behind the stylized facts is not fully 
understood. In order to imitate the real financial markets and to understand the 
origins of the stylized facts, a variety of models have been proposed a nd e xamined. 
Actually many models are able to capture some of the stylized factPH. 

In empirical finance the volatility is an important value to measure the risk. One 
of the stylized facts of the volatility is that the volatility of price returns changes 
in time and shows clustering, so called "volatility clustering". Then the histogram 
of the resulting price returns shows a fat-tailed distribution which indicates that 
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the probabihty of having a large price change is higher than that of the Gaussian 
distribution. In order to mimic these empirical properties of the volatility and to 
forecast the future volatility values, Engle advocated the autoregressive conditional 
hetroskedasticity (ARCH) model^ where the volatility variable changes determin- 
istically depending on the past squared value of the return. Later the ARCH model 
is generalized by adding also the past volatility dependence to the volatility change. 
This model is known as the generalized ARCH (GARCH) modeP^. The parameters 
of the GARCH model applied to financial time series are conventionally determined 
by the maximum likelihood method. There are m any ex tended versions of GARCH 
models, such as EGARClil^^ GJ'^, QGARClil912Q! models etc., which are de- 
signed to increase the ability to forecas t the v olatility value. 

The stochastic volatility (SV) modern^ is another model which captures the 
properties of the volatility. In contrast to the GARCH model, the volatility of the SV 
model changes stochastically in time. As a result the likelihood function of the SV 
model is given as a multiple integral of the volatility variables. Such an integral in 
general is not analytically calculable and thus the determination of the parameters 
of the SV model by the maximum likelihood method becomes difficult. To overcome 
this difficulty in the maximum likelihood method the Markov Chain Monte Carlo 
(MCMC) method based on the Bayesian approach is proposed and develope cPI. In 
the MCMC of the SV model one has to update not only the parameter variables 
but also the volatility ones from a joint probability distribution of the parameters 
and the volatility variables. The number of the volatility variables to be updated 
increases with the data size of time series. The first proposed update scheme of 
the volatility variables is based on the local update such as the Metropolis-type 
algorithni^-'^'. It is however known that when the local update scheme is used for 
the volatility variables having interactions to their neighbor variables in time, the 
autocorrelation time of sampled volatility variables becomes large and thus the local 
update scheme becomes ineffective In order to improve the efficiency of the local 
update method the blocked scheme which updates several variables at once is also 
proposecPSEH, A recent survey on the MCMC studies of the SV model is seen in 
Ref.25. 

In our study we use the HMC algorithnl^ which had not been considered 
seriously for the MCMC simulation of the SV model. In ffnance there exists an 
application of the HMC algorithm to the GARCH modeP3 where three GARCH 
parameters are updated by the HMC scheme. It is more interesting to apply the 
HMC for updates of the volatility variables because the HMC algorithm is a global 
update scheme which can update all variables at once. This feature of the HMC 
algorithm can be used for the global update of the volatility variables whi ch c an not 
be achieved by the standard Metropolis algorithm. A preliminary studjPSl shows 
that the HMC algorithm samples the volatility variables effectively. In this paper we 
give a detailed description of the HMC algorithm and examine the HMC algorithm 
with artiffcial financial data up to the data size of T=5000. We also make an 
empirical analysis of the Nikkei 225 stock index by the HMC algorithm. 
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2. Stochastic Volatility Model 

The standard version of the SV modePll^ is given by 

Ut = ctet = exp(/it/2)et, (1) 

ht = iJ. + (i){ht-i- +rit, (2) 

where yt = {yi, y2, ■■■,yn) represents the time series data, ht is defined by ht = Ina^ 
and at is called volatility. We also call ht volatility variable. The error terms et and rjt 
are taken from independent normal distributions -/V(0, 1) and ^'^(O, cr^) respectively. 
We assume that |<^| < 1. When (f) is close to one, the model exhibits the strong 
persistency of the volatility shock. 

For this model the parameters to be determined are /i, (j) and a^. Let us use 9 
as 9 — (/i, (/), (7^). Then the likelihood function L{6) for the SV model is written as 

/n 
n f{et\a^)fiht\9)dhidh2...dhn, (3) 

where 

fiet\af) = {27ra^)-Kxp(-^y (4) 

/(/..|D) = (2,,;)-'.xp(- ""-''-^^^--'-'l' ). (6) 

As seen in Eq.Q, L{9) is constructed as a multiple integral of the volatility vari- 
ables. For such an integral it is difficult to apply the maximum likelihood method 
which estimates values of 9 by maximizing the likelihood function. Instead of using 
the maximum likelihood method we perform the MCMC simulations based on the 
Bayesian inference as explained in the next section. 



3. Bayesian inference for the SV model 

From the Bayes' rule, the probability distribution of the parameters 9 is given by 

f{9\y) ^ ^L{9)n{9), (7) 

where Z is the normalization constant Z = j L{9)'K{9)d9 and Tr{9) is a prior disti- 
bution of 9 for which we make a certian assumption. The values of the parameters 
arc inferred as the expectation values of 9 given by 

(9) - J 9f{9\y)d9. (8) 

In general this integral can not be performed analytically. For that case, one can 
use the MCMC method to estimate the expectation values numerically. 
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In the MCMC method, we first generate a series of 9 with a probability of 
P{9) ^ f{e\y). Let 6I« = (6l(i), 6'(2), 6'(''')) be values of 6* generated by the MCMC 
sampling. Then using these k values the expectation value of 9 is estimated by an 
average as 

k 

(^?> = i^0«. (9) 

4=1 

The statistical error for k independent samples is proportional to —=. When the 

. [2^ 

sampled data are correlated the statistical error will be proportional to \ / — where 

V k 

T is the autocorrelation time between the sampled data. The value of r depends 
on the MCMC sampling scheme we take. In order to reduce the statistical error 
within limited sampled data it is better to choose an MCMC method which is able 
to generate data with a small t. 



3.1. MCMC Sampling of 9 

For the SV model, in addition to 9, volatility variables ht also have to be updated 
since they should be integrated out as in Eq.Q. Let P{9, ht) be the joint probability 
distribution of 9 and ht- Then P{9, ht) is given by 

Pi9,ht)^ Li9,ht)n{9), (10) 

where 

n 

L{9,ht)^Y[f{et\ht)f{ht\9). (11) 
t=i 

For the prior tt{9) we assume that 7r(cr^) ~ (c^)^^ and for others tt{h) = TT{(f>) ~ 
constant. ^^^^ 

The MCMC sampling methods for 9 are given in the followinJ2Il22]^ r^j^^ prob- 
ability distribution for each parameter can be derived from Eq. (jlO|) by extracting 
the part including the corresponding parameter. 

• (T^ update scheme. 

The probability distribution of is given by 

^('x^)-(^^)-^-'exp(^-A^, (12) 

where 

n 

A = ^{{l - 4>^){h, - fi)' + J2[ht - M - ^{ht-i - ^^)?}■ (13) 

t=2 

Since Ea. p2p is an inverse gamma distribution we can easily draw a value 
of a^- by using an appropriate statistical library in the computer. 
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• fi update scheme. 

The probabihty distribution of /i is given by 

^(M)-exp|-^0i-|)2|, (14) 

where 

B={l-^^) + {n-l){l~cl,r, (15) 

and 

n 

C = (1 - cl)^)hi + (1 - 0) Y^iht - 4>ht-i)- (16) 

t=2 

fi is drawn from a Gaussian distribution of Eq. (I14p . 

• (j) update scheme. 

The probabihty distribution of cj) is given by 

P(</)) ^ (1 - 02)1/2 ,^p|__^(^ _ |)2^^ (^7) 

where 

n 

D = -(/ii - + ^(/it_i - /z)2, and £; = EtiC^t - M)(/it-i - m). (18) 



In order to u pdate cj) 'with Eq. pTj) . we use the Metropohs-Hastings 
algorithnJSnED. Let us write Eq.([T71) as P{(t)) ^ Pi{4>)P2{4>) where 

PM)^{\-^^Y'\ (19) 
P2(0)^exp{-^(0-|)2}. (20) 

Since P2{<t>) is a Gaussian distribution we can easily draw from Eq. (|20p . 
Let (pnew be a candidate given from Ea. (|20p . Then in order to obtain the 
correct distribution, ipnew is accepted with the following probability Pmh- 



Pmh = min 



In addition to the above step we restrict cj) within [—1, 1] to avoid a negative 
value in the calculation of square root. 

3.2. Probability distribution for ht 

The probability distribution of the volatility variables ht is given by 

P{ht)=P{huh2,...,K)^ (22) 

y l^i=l\ 2 ~^ 2^ J 2ay(l-4,-i) l^i=2 2^ ) ' 
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This probability distribution is not a sim ple fun ction to draw values oi ht- A conven- 
tional method is the Metr opolis meth oc j'^^ l '^^ l which updates the variables locally. 
There are several developed to update ht from Eq. ((22| . Here 

we use the HMC algorithm to update ht globally. The HMC algorithm is described 
in the next section. 



4. Hybrid Monte Carlo Algorithm 

Originally the HMC algorithm is developed for the MCMC simulations of the lattice 
Quantum Chromo Dynamics (QCD) calculations^'^. A major difhculty of the lattice 
QCD calculations is the inclusion of dynamical fermions. The effect of the dynamical 
fermions is incorporated by the determinant of the fermion matrix. The computa- 
tional work of the determinant calculation requires 0(y^) arithmetic operationi^Slj 
where V is the volume of a 4-dimensional lattice. A typical size of the volume is 
V > 10"'. The standard Metropolis algorithm which locally updates variables does 
not work since each local update requires 0{V'^) arithmetic operations for a deter- 
minant calculation, which results in unacceptable computational cost in total. Since 
the HMC algorithm is a global update method, the computational cost remains in 
the acceptable region. 

The basic idea of the HMC algorithm is a combination of molecular dynamics 
(MD) simulation and Metropolis accept /reject step. Let us consider to evaluate the 
following expectation value {0{x)) by the HMC algorithm. 

(0(x)> / 0{x)f{x)dx = / Oix)e^"f^'^'>dx, (23) 



where x — {xi,X2, ■■■,Xn), f(x) is a probability density and 0{x) stands for an 
function of x. First we introduce momentum variables p — (pi,P2, ■■■,Pn) conjugate 
to the variables x and then rewrite Eq. (l23|) as 

{0{x)) = ^ J 0{x)e-2P'+i''fi^)dxdp=^ J 0{x)e-"^P'-'Uxdp. (24) 

where Z is a normalization constant given by 

Z= fcxp (-IpA dp, (25) 



2 

and p'^ stands for pf- H{p, x) is the Hamiltonian defined by 

H{p,x)^^p^ -Infix). (26) 

Note that the introduction oi p docs not change the value of {0{x)). 

In the HMC algorithm, new candidates of the variables {p, x) are drawn by 
integrating the Hamilton's equations of motion, 

dxi _ dH 
dpi _ dH 
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In general the Hamilton's equations of motion arc not solved analytically. Therefore 
we solve them numerically by doing the MD simulation. Let TMoi^t) be an elemen- 
tary MD step with a step size At, which evolves {p{t),x{t)) to + At), a;(t + At)): 



-MD 



(At) : (p(t), x{t)) -> {pit + At),x{t + At)). (29) 



Any integrator can be used for the MD simulation provided that the following 
conditions are satisfiecPSl 

• area preserving 

dp{t)dx{t)dx = dp{t + At)dx{t + At). (30) 

• time reversibility 

TMoi-At) : {p{t + At), x{t + At)) ^ (p(t), x{t)). (31) 

The simplest and often used integrator satisfying the above two conditions is 
the 2nd order leapfrog integrator given by 

x,(t + At/2) = x,{t) + ^p^(t) 

K(t + At) =p(t), -At— 

x^(t + At) = .T4t + At/2) + ^p,;(t + At). (32) 

In this study we use this integrator. The numerical integration is performed TV steps 
repeatedly by Eg. ((5^ and in this case the total trajectory length A of the MD is 
\ = N X At. 

At the end of the trajectory we obtain new candidates {p' , x'). These candidates 
are accepted with the Metropohs test, i.e. {p',x') are globally accepted with the 
following probability, 

where AH is the energy difference given by AH — H{p', x') — H{p, x). Since we inte- 
grate the Hamilton's equations of motion approximately by an integrator, the total 
Hamiltonian is not conserved, i.e. AH ^ 0. The acceptance or the magnitude of AH 
is tuned by the step size At to obtain a reasonable acceptance. Actually th ere ex- 
ists the optimal acceptance which is about 60 — 70% for 2nd order integrator J32l33]_ 
Surprisingly the optimal acceptance is not dependent of the model we consider. For 
the n-th order integrator the optimal acceptance is expected to b J^S! ^ gxp ( 

We could also use higher order integrators which give us a smaller energy dif- 
ference AH. However the higher order integrators are not always eff ective since 
they need more arithmetic operations than the lower order integratorJSHSSl xhe 
efficiency of the higher order integrators depends on the model we consider. There 
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also exist improved integrators which have less arithmetic operations than the con- 
ventional integratorj2ll. 

For the volatility variables ht, from Eq. p2|) . the Hamiltonian can be defined by 

HiPu h) - + 1+2^27(13^+2. 2^2 ' (34) 

i=l 4=1 ^' ^ ^ ' i=2 ^ 

where pi is defined as a conjugate momentum to hi. Using this Hamiltonian we 
perform the HMC algorithm for updates of ht . 



5. Numerical Studies 

In order to test the HMC algorithm we use artificial financial time series data 
generated by the SV model with a set of known parameters and perform the MCMC 
simulations to the artificial financial data by the HMC algorithm. We also perform 
the MCMC simulations by the Metropolis algorithm to the same artificial data and 
compare the results with those from the HMC algorithm. 

Using Eq.([T]) with (/> — 0.97,(t^ — 0.05 and fi — —1 we have generated 5000 
time series data. The time series generated by Eq.([T]) is shown in FiglT] From those 
data we prepared 3 data sets: (1)T=1000 data (the first 1000 of the time series), 
(2)T=2000 data (the first 2000 of the time series) and (3) T=5000 (the whole data). 
To these data sets we made the Bayesian inference by the HMC and Metropolis 
algorithms. Precisely speaking both algorithms are used only for the MCMC update 
of the volatility variables. For the update of the SV parameters we used the update 
schemes in Sec.3.1. 

For the volatility update in the Metropolis algorithm, we draw a new candidate 
of the volatility variables randomly, i.e. a new volatility /i"*^"" is given from the 
previous value hf"^ by 

h^<^^ = hf"^ + 6{r~-0.5), (35) 

where r is a uniform random number in [0, 1) and (5 is a parameter to tune the 
acceptance. The new volatility /i"'^™ is accepted with the acceptance Pmetro 



new 



(36) 



p{hf'^) 

where P{ht) is given by Eq. (|22|) . 

The initial parameters for the MCMC simulations are set to (f> — 0.5,cr,'j = 1.0 
and = 0. The first 10000 samples are discarded as thermalization or burn- in 
process. Then 200000 samples are recorded for analysis. The total trajectory length 
A of the HMC algorithm is set to A = 1 and the step size Ai is tuned so that the 
acceptance of the volatility variables becomes more than 50%. 

First we analyze the sampled volatility variables. Figl^] shows the Monte Carlo 
(MC) history of the volatility variable /iioo from T — 2000 data set. We take /iioo 
as the representative one of the volatility variables since we have observed the 
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t 



Fig. 1. The artificial SV time series used for tliis study. 




55000 

Monte Carlo history 




55000 

Monte Carlo history 



Fig. 2. Monte Carlo histories of hioo generated by HMC (left) and Metropolis (right) with 
T = 2000 data set. The Monte Carlo histories in the window from 50000 to 60000 are shown. 



similar behavior for other volatility variables. See also Fig[3]for the similarity of the 
autocorrelation functions of the volatility variables. 

A comparison of the volatility histories in Fig[2]clearly indicates that the correla- 
tion of the volatility variable sampled from the HMC algorithm is smaller than that 
from the Metropolis algorithm. To quantify this we calculate the autocorrelation 
function (ACF) of the volatility variable. The ACF is defined as 

wEf-Mj) - {x)){x{j + t) - {x}) 

AC Fit) = ^^'-'^ ^ ' -4^—^ (37) 

where {x) and cr^ are the average value and the variance of x respectively. 

Fig [3] shows the ACF for three volatility variables, ft.io,/i20 and ft-ioo sampled 
by the HMC. It is seen that those volatility variables have the similar correlation 
behavior. Other volatility variables also show the similar behavior. Thus hereafter 
we only focus on the volatility variable /iioo as the representative one. 

FigHl compares the ACF of hioo by the HMC and Metropolis algorithms. It 
is obvious that the ACF by the HMC decreases more rapidly than that by the 
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20 40 60 80 

t 

Fig. 3. Autocorrelation functions of three volatility variables /iio,/i20 a-^d /iioo sampled by the 
HMC algorithm for T = 2000 data set. These autocorrelation functions show the similar behavior. 




100 200 300 400 500 

t 

Fig. 4. Autocorrelation function of the volatility variable /iioo by the HMC and Metropolis 
algorithms for T = 2000 data set. 

Metropolis algorithm. We also calculate the autocorrelation time Tint defined by 

1 °° 

nnt^ ^+Y.ACFit). (38) 

The results of Tint of the volatility variables are given in Table 1. The values in 
the parentheses represent the statistical errors estimated by the jackknife method. 
We find that the HMC algorithm gives a smaller autocorrelation time than the 
Metropolis algorithm, which means that the HMC algorithm samples the volatility 
variables more effectively than the Metropolis algorithm. 

Next we analyze the sampled SV parameters. Figl5] shows MC histories of the 
parameter sampled by the HMC and Metropolis algorithms. It seems that both 
algorithms have the similar correlation for (/>. This similarity is also seen in the ACF 
in Figint^left), i.e. both autocorrelation functions decrease in the similar rate with 
time t. The autocorrelation times of (p are very large as seen in Table 1. We also 
find the similar behavior for cr^, i.e. both autocorrelation times of are large. 

On the other hand we see small autocorrelations for fi as seen in FiglSf right). 
Furthermore we observe that the HMC algorithm gives a smaller Tint for fi than 
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M 




hioo 


true 


0.97 


-1 


0.05 




T=1000 


HMC 


0.973 


-1.13 


0.053 






SD 


0.010 


0.51 


0.017 






SE 


0.0004 


0.003 


0.001 








360(80) 


3.1(5) 


820(200) 


12(1) 




Metropolis 


0.973 


-1.14 


0.053 






SD 


0.011 


0.40 


0.017 






SE 


0.0005 


0.003 


0.0013 








320(60) 


10.1(8) 


720(160) 


190(20) 


T=2000 


HMC 


0.978 


-0.92 


0.053 






SD 


0.007 


0.26 


0.012 






SE 


0.0003 


0.001 


0.0009 








540(60) 


3(1) 


1200(150) 


18(1) 




Metropolis 


0.978 


-0.92 


0.052 






SD 


0.007 


0.26 


0.011 






SE 


0.0003 


0.003 


0.0009 






2Ti„t 


400(100) 


13(2) 


1000(270) 


210(50) 


T=5000 


HMC 


0.969 


-1.00 


0.056 






SD 


0.005 


0.11 


0.009 






SE 


0.0003 


0.0004 


0.0007 








670(100) 


4.2(7) 


1250(170) 


10(1) 




Metropolis 


0.970 


-1.00 


0.054 






SD 


0.005 


0.12 


0.008 






SE 


0.00023 


0.0011 


0.0005 






2ri„t 


510(90) 


30(10) 


960(180) 


230(28) 



Table 1. Results estimated by the HMC and Metropolis algorithms. SD stands for Standard 
Deviation and SE stands for Statistical Error. The statistical errors are estimated by the jackknife 
method. We observe no significant differences on the autocorrelation times among three data sets. 



that of the Metropolis algorithm, which means that HMC algorithm samples fi 
more effectively than the Metropolis algorithm although the values of Tint for jj, 
take already very small even for the Metropolis algorithm. 

The values of the SV parameters estimated by the HMC and the Metropolis 
algorithms are listed in Table 1. The results from both algorithms well reproduce 
the true values used for the generation of the artificial financial data. Furthermore 
for each parameter and each data set, the estimated parameters by the HMC and 
the Metropolis algorithms agree well. And their standard deviations also agree 
well. This is not surprising because the same artificial financial data, thus the same 
likelihood function is used for both MCMC simulations by the HMC and Metropolis 
algorithms. Therefore they should agree each other. 
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45000 
MC history 



45000 
MC history 



50000 



Fig. 5. Monte Carlo histories of </> generated by HMC (left) and Metropolis (right) for T = 2000 
data set. 




300 



Fig. 6. Autocorrelation functions of (/((left) and /^(right) by the HMC and Metropolis algorithm 
for T = 2000 data set. 



6. Empirical Analysis 

In this section we make an empirical study of the SV model by the HMC algorithm. 
The empirical study is based on daily data of the Nikkei 225 stock index. The 
sampling period is 4 January 1995 to 30 December 2005 and the number of the 
observations is 2706. FiglTfleft) shows the time series of the data. Let pi be the 
Nikkei 225 index at time i. The Nikkei 225 index pi are transformed to returns as 



= 1001n(p,/p, 



(39) 



where s is the average value of hi{pi/pi^i). FiglTl^right) shows the time series of 
returns calculated by Eq.l[39|). We perform the same MCMC sampling by the HMC 
algorithm as in the previous section. The first 10000 MC samples are discarded and 
then 20000 samples are recorded for the analysis. The ACF of sampled /iioo and 
sampled parameters are shown in Fig[5] Qualitatively the results of the ACF are 
similar to those from the artificial financial data, i.e. the ACF of the volatility and 
/i decrease quickly although the ACF of (p and cr^ decrease slowly. The estimated 
values of the parameters are summarized in Table 2. The value of (j) is estimated to 
he (p 0.977. This value is very close to one, which means the time series has the 
strong persistency of the volatility shock. The similar values are also seen in the 
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HMC 


4^ 










0.977 


0.52 


0.020 




SD 


0.006 


0.13 


0.005 




SE 


0.001 


0.0016 


0.001 






560(190) 


4(1) 


1120(360) 


21(5) 



25000 



20000 



15000 



10000 



Table 2. Results estimated by the HMC for the Nikkei 225 index data. 



Nikkei 225 Index 




u" 




500 1000 1500 2000 2500 3000 
t 



500 1000 1500 2000 2500 3000 
t 



Fig. 7. Nikkei 225 stock index from 4 January 1995 to 30 December 2005{left) and returns (right). 




200 400 600 800 1000 
t 



Fig. 8. Autocorrelation functions of the volatility variable /iioo (left) and the sampled parameters 
(right). 



previous 



7. Conclusions 

We applied the HMC algorithm to the Bayesian inference of the SV model and 
examined the property of the HMC algorithm in terms of the autocorrelation times 
of the sampled data. We observed that the autocorrelation times of the volatility 
variables and /i parameter are small. On the other hand large autocorrelation times 
are observed for the sampled data of (j) and af^ parameters. The similar behavior 



for the autocorrelation times are also seen in the literature 

From comparison of the HMC and Metropolis algorithms we find that the HMC 
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algorithm samples the volatility variables and fi more effectively than the Metropolis 
algorithm. However there is no significant difference for and sampling. Since 
the autocorrelation times of /i for both algorithms are estimated to be rather small 
the improvement of sampling /i by the HMC algorithm is limited. Therefore the 
overall efficiency is considered to be similar to that of the Metropolis algorithm. 

By using the artificial financial data we confirmed that the HMC algorithm cor- 
rectly reproduces the true parameter values used to generate the artificial financial 
data. Thus it is concluded that the HMC algorithm can be used as an alternative 
algorithm for the Bayesian inference of the SV model. 

If we are only interested in parameter estimations of the SV model, the HMC 
algorithm may not be a superior algorithm. However the HMC algorithm samples 
the volatility variables effectively. Thus the HMC algorithm may serve as an efficient 
algorithm for calculating a certain quantity including the volatility variables. 
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